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, We evaluate the angular bispectrum of the CMB temperature anisotropy at large angular scale due 

^SJ ' to a stochastic background of primordial magnetic fields. The shape of non-Gaussianity depends on 

the spectral index of the magnetic field power spectrum and is peaked in the squeezed configuration 
for a scale-invariant magnetic spectrum. By using the large angular part of the bispectrum generated 
by magnetic fields, the present bounds on non-Gaussianity set a limit on the amplitude of the 
primordial magnetic field of the order of 0(10) nGauss for the scale-invariant case and 0(20) nGauss 
for the other spectral indexes. 



in 



I 

o 



O, 

O . I. INTRODUCTION 

Cl^' Cosmological inflation [l| has become the dominant paradigm to understand the initial conditions for the Cosmic 
Microwave Background (CMB) anisotropies and structure formation. This picture has recently received further 
spectacular confirmation by the Wilkinson Microwave Anisotropy Probe (WMAP) five year set of data Present 
^ ' [^] and future [3| experiments may be sensitive to the non-linearities of the cosmological perturbations at the level 
d , of second- or higher-order perturbation theory. The detection of these non-linearities through the non-Gaussianity 
(NG) in the CMB 3 

has become one of the primary experimental targets. 
, A possible source of NG could be primordial in origin, being specific to a particular mechanism for the generation 
J> ■ of the cosmological perturbations. This is what makes a positive detection of NG so relevant: it might help in 
discriminating among competing scenarios which otherwise might be indistinguishable. Indeed, various models of 
04 inflation, firmly rooted in modern particle physics theory, predict a significant amount of primordial NG generated 
\l either during or immediately after inflation when the comoving curvature perturbation becomes constant on super- 
" I horizon scales [^]. While single-field Q and two(multi)-field [6] models of inflation generically predict a tiny level of 
NG, 'curvaton-type models', in which a significant contribution to the curvature perturbation is generated after the 
, end of slow-roll inflation by the perturbation in a field which has a negligible effect on inflation, may predict a high 
■ level of NG 0, [1| • Alternatives to the curvaton model are those models characterised by the curvature perturbation 
being generated by an inhomogeneity in the decay rate [l3| or the mass [ll| or of the particles responsible for 
the reheating after inflation. Other opportunities for generating the curvature perturbation occur at the end of 
inflation [12] and during preheating JJ|J. All these models generate a level of NG which is local as the NG part of the 
primordial curvature perturbation is a local function of the Gaussian part, being generated on superhorizon scales. 
' In momentum space, the three point function, or bispectrum, arising from the local NG is dominated by the so-called 
"squeezed" configuration, where one of the momenta is much smaller than the other two and it is parametrized by the 
non-linearity parameter /^l- Other models, such as DBI inflation [l3| and ghost inflation Jl^, predict a different kind 
of primordial NG, called "equilateral", because the three-point function for this kind of NG is peaked on equilateral 
configurations, in which the lengths of the three wave- vectors forming a triangle in Fourier space are equal [l^. The 
equilateral NG is parametrized by an amplitude f^^^ [TtS- Present limits on NG are summarised by —9 < /^l < HI 
and -151 < /nI;'" < 253 at 95% CL 
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On the other hand there might exist other sources of primordial NG in the CMB anisotropies beyond the primordial 
ones related to the dynamics of inflation. One interesting possibility is the contribution to the non-Gaussian signal 
in the CMB anisotropies from a stochastic background of primordial magnetic fields. Large scale magnetic fields 
are almost everywhere in the universe, from galaxies up to those present in galaxy clusters and in the inter-cluster 
medium [19 ]. The dynamo effect provides a mechanism to explain the observed magnetic fields associated to galaxies, 
whereas those associated to clusters may be generated by gravitational compression. Both these mechanisms require 
an initial magnetic seed, although with different amplitude and different correlation length. 

Possible explanations for this initial seed have driven the interest in primordial magnetic fields generated in the 
early universe. A stochastic background of primordial magnetic fields (PMFs) generated in the early universe with 
a mean amplitude well below micro-Gauss level can leave imprints on the temperature and polarisation anisotropy 
pattern of the cosmic microwave background (CMB). The impact of a stochastic background of PMFs onto CMB 
anisotropies has distinctive imprints, such as a contribution in temperature which is larger than the CMB angular 
power spectrum sourced by scalar cosmological perturbation at high i and a contribution in polarisation which include 
either BB (generated by vector and tensor perturbations or by Faraday rotation [l^l) or parity-odd correlators as TB 
(generated by an helical component [2ll|). 

As we mentioned, another distinctive imprint of PMF in CMB anisotropies is its non-Gaussian nature. The CMB 
signature of this type which has been first considered in the literature is due to a homogeneous PMF. A homogeneous 
magnetic field with fixed direction breaks spatial isotrop y in the universe, and therefore leads to non-zero correlations 
between multipoles at different {a'i-i,m0.i+i,m) 7^ [22j |. This effect has been first proposed in [2^, and arises 
through the generation of vector metric pertubations from the Alfven waves magnetically induced in the primordial 
fluid. Recently, it has been reanalysed and found to reproduce, for a sufficiently high magnetic field amplitude, some 
of the anomalies of the CMB large scale fluctuations observed by WMAP such as the north-south asymmetry and the 
quadrupole-octopole alignment [24| . 

This effect is related to the presence of a homogeneous magnetic field (or equivalently, a stochastic magnetic field 
with correlation length larger than the horizon today). On the other hand, the CMB contribution of a stochastic 
background of PMFs, modelled as a fully inhomogeneous component, is intrinsically non-Gaussian: the PMFs energy- 
momentum tensor, the Lorentz force acting on baryons are quadratic in the magnetic field B(?7, x), which is randomly 
distributed with a Gaussian distribution function. The source terms to the Einstein-Boltzmann system are therefore 
X-distributed, leading to a PMF contribution to CMB fully non-Gaussian. Higher order statistical moments of the 
energy-momentum tensor of PMFs are therefore non-vanishing at leading order and are calculable with minimal 
assumptions, such as cutting sharply the power spectrum beyond a certain scale ko [25j . 

Non-gaussianities from PMFs are much less studied than those generated in inflationary methods. The study of 
the three point statistics of the PMF energy-momentum tensor in [25| is limited to the simplest particular coUincar 
configuration. Nevertheless, due to the presence of a contribution from the coUinear configuration, one deduces 
that non-gaussianities from PMFs can be different from the inflationary case, in which the coUinear contribution is 
generically negligible with respect to the equilateral and squeezed ones (see however [1^). In this paper we focus on 
the three point statistics of the PMF energy density, studying the contribution of all three configurations. By using 
the large scale relation between temperature anisotropies and PMF energy density given in Ref. 21] , we compute the 
temperature bispectrum and compare its contribution to the non-Gaussian statistics in the CMB anisotropies with 
the present observational bounds. We also compare our results with those of the very recent paper [28l|. 

Our paper is organised as follows. In section II we introduce the stochastic background of primordial magnetic 
fields and discuss the infrared behaviour of the spectra of its energy- density. In section III and IV we discuss the 
CMB temperature spectrum and bispectrum induced by PMF on large scales. Section V is devoted to the analytic 
computation of the magnetic energy density bispectrum (/9B(k)p_B(q)pB(p)) in general and for the coUinear, squeezed 
and equilateral configurations. In Section VI we insert these results into the CMB temperature bispectrum on large 
scales, and we give an estimation of the signal in section VII. In the first Appendix we derive analytic approximations 
to some integrals of Bessel functions which are useful to calculate both the spectrum and the bispectrum, and in the 
second Appendix we give the details for the exact computation of the energy density bispectrum in the coUinear case 
for n = 2, -2. 



II. PRIMORDIAL STOCHASTIC MAGNETIC FIELD 



We adopt notations consistent with [2^, [s^l : 
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where the definition of the delta function comes from J (fikj {2ttY e *'''^(5(k) = l/(27r)'^. With these conventions, the 
magnetic field power spectrum (defined as the Fourier transform of the two point correlation function) is^: 

(B.(k)B*(q)) = {2Tif5\\<i-ci){5,j~hk)PB{k) (1) 
Ps(fc) ^ A¥\ k<kD, (2) 

where ki = ki/k, A is a normalisation constant, n the spectral index and kjj the upper cutoff. Using the above 
equations we can define the mean square of the magnetic field as 

(i?2(x)) = 4^. (3) 

If we are interested in the mean amplitude of the magnetic field on a given characteristic scale, we smooth the power 
spectrum over the chosen scale using a Gaussian filter: we have then {B'^{x))\\ — B\ with 



so that 



2 n + 'i ^,n + i. 



We also define the adimensional quantity il^' given by the ratio of the magnetic and the total radiation energy 
densities: 

"S^Prci" (10-9Gauss)2 ' 

where for the last equality we have used Prci(?7o) — 2 x 10~^^ GeV'', and 770 denotes the conformal time today. 

The upper cutoff fc^i corresponds to the dam ping sc ale, representing the dissipation of magnetic energy due to 
the generation of magneto-hydrodynamic waves [3ll |32| . Alfven waves are the most effective in dissipating magnetic 
energy, and in [3^ it is demonstrated that around recombination the damping occurs at scales k~^ <, k~^ ~ VAisiik, 
where Va is the Alfven speed and Lsiik the Silk damping scale at recombination. Strictly speaking, Alfven waves 
are oscillatory perturbations superimposed on a homogeneous magnetic component, and the Alfven speed depends 
on the amplitude of the homogeneous component. In the cosmological context where the magnetic field is purely 
stochastic, the amplitude of this component can be taken as the one of a 'low frequency' component obtained by 
smoothing the magnetic field amplitude over a sufficiently large scale [s^l- This scale corresponds to the Alfven 
scale at recombination, k~^ ~ VAr]rec' magnetic modes on lager scales, in fact, do not have time to oscillate before 
recombination jsl]- One has therefore k^/kA = ?7rec/isiik — 10, and = Bj^^/{4:n{p + p)). Consequently, the upper 
cutoff at the epoch of recombination is given by 



° LsiikY 3 Bl^ Lsiik Vf^V (" + 3)r((" + 3)/2) ' 

Finiteness of the total magnetic energy density implies n > — 3. In the rest of the paper we keep the spectral index 
as a free parameter, when possible; however, in order to carry on our calculations analytically we sometimes need to 
specify it. For example, in Section |Vl we choose the values n = 2 and n = —2. n = 2 is the lowest possible value 
for a magnetic field generated by a causal process [s^ , such as a phase transitio n [35I , or a charge separation process 
operating around recombination [36|| . A magnetic field generated during inflation |37l| (or by any other a-causal process 
such as, for example, in pre big bang theories [38|]), can take any value of n > —3. However, because of Nucleosynthesis 
constraints [1^ , only for very red spectra the magnetic field can assume sufficiently high amplitudes to have an impact 
in the CMB. Therefore, in the following we choose the value n = —2 (for some analytic calculations), or n —3 
when possible. 



In this paper we neglect the possible presence of an helical component for the magnetic field, see for example (2 ill 
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The spatial part of the energy momentum tensor of the magnetic field is 



1 



47r 



and the magnetic energy density 



i(5,yB,(x)B,(x)-B,(x)Bj(x) 



1 



PS(X) = T,,(X) = ^B2(x), 
P_ 

{2^f 



P-^"^^ ^ 8l/(S^^^('^-P)^^(P)• 



(8) 

(9) 
(10) 



As we will see in the next section, to calculate the CMB temperature spectrum from the Sachs Wolfe effect we need 
the power spectrum of the magnetic energy density: 



(pB(kK(q))^(27r)35(k-q)K(fc)|^ 



(8^)^ 



,,5(k-q) j d3pPB(p)PB(|k-p|)(l + /i'), 



the second equality is obtained using Eq. PTI)1 . and /i = p • k — p. Therefore 



\PB{kf 



/ rf'P^B(p)PB(|k-p|)(l + M'). 



(11) 



(12) 



As demonstrated in Ref. [27l|, |pB(fc)P goes to zero at a wavenumber corresponding to twice the magnetic field 
spectrum cutoff, k — 2kD- Eq. (21) of [27| gives the behaviour of \PB{k)\^ at large scales k <^ kn and for spectral 
indexes n > —3/2: the generic behaviour in this case is white noise [27| 



\PBik)\' 



327r4(2n + 3) ' 

For n = —2, an exact calculation as in Refs. [13,14^ gives the behaviour for |pb(/c)P at large scales ^ fc/j as 

"^^(^)"'"5T2^' 
For n < — 3/2 we use the approximated formula given by Ref. [29| : 

3A^ n 



\PB{k)\^ 



L2n+3 



1287r4 (2n + 3)(7i + 3) 



(13) 



(14) 



(15) 



For n = —2, expressions (|14p and (jlSp are in agreement concerning the dependence on the wavenumber, but the 
numerical factor differs by a factor tt^/S. 



III. CMB TEMPERATURE SPECTRUM AT LARGE ANGULAR SCALES 

We use the characterisation of the CMB temperature anisotropy induced by a stochastic background of primordial 
magnetic fields given in It is useful to define the adimensional magnetic energy parameter in fc-space 

0^(k) = ^. (16) 

From the initial conditions given in [2?] . [40| we parametrize the temperature anisotropy as 

ls, + ^p^^nBik), (17) 

where a ~ 0.1 is a multiplication constant required since the above equation would be exact with a = 1 in the 
radiation era. Therefore, the temperature anisotropy is given in terms of this quantity as 

"y^"; ' = ^ilsikMkilJo - %ec)) , (18) 
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where ji is the spherical Bessel function and 770, ?7doc denote conformal time respectively today and at decoupling. 
The CMB power spectrum is therefore [30]: 



TT Jo (2^+1)2 Svr 

In the case n > —3/2, substituting definition and Eq. in the above equation we have [ij 



8 2^52n + 3p2 fc^?7o^io 5127r 2n + 3 fc^^*^ 



where a; = fcryo, we have approximated ji{k{rjQ — jydoc)) — jeikrio) and we integrate only up to the upper cutoff kjj since 
we are using the approximated expression Eq. (|I3p which is strictly valid only for k <^ k^,. For the second equality 
in the above equation, we have approximated the integral as given in Eq. (jA2|) of appendix [Aj since we have that 
y — korjo » 1. We have also used Eq. ^ to express the result in terms of the mean squared magnetic field. 
For n = —2, we use instead Eq. (fT4|) : substituting it in Eq. p9|) . we find 



W - ,nn^ ^ 2 — ^ dxxjf{x) ~ ^ ^og [ ^— ] for n = -2 , (21) 



men^pl^ikDVo)^ Jo '~ 8192 p^, (fc^^) 

where in the second equality we use the approximation given in Eq. (|A4[) of appendix |^ For more negative values of 
n, n < —2, in the absence of an exact expression, we use Eq. (|15p : the CMB spectrum becomes 

3 0-2 /12 l,2n+6 /-fcoiJo 

- 1024 7r5p2^i(fcz5?7o)2"+M2n + 3)(n + 3)yo ^ ^ 

~. 3.2 „(. + 3) r[-n-2] (i.2^^ 1 ^^^^^^ ^^^^^^^^ ^^^^ 



4096V^(2n + 3)r[-n-3/2] p^, (fcc77o)2"+6 

where again for the second equality we have used Eq. (jA3| of appendix [A] The slope in £ of this last expression is 
only approximatively recovered numerically, whereas it is perfectly recovered for Eqs. (|20l2ip . 



IV. CMB TEMPERATURE BISPECTRUM AT LARGE ANGULAR SCALES 

We want to evaluate the CMB angular bispectrum of the temperature anisotropy due to the Sachs Wolfe effect 
induced by the magnetic field energy density. The angular bispectrum is given by {ai-^rni(ie2m2(ii3m3) , with the 
spherical harmonic expansion coefficients 

a,„(x)= [ dnnY;^{h;e)e^°\^,h), (23) 



where Yg*^(h;e) is the spherical harmonic with respect to a basis where e is an arbitrary but fixed direction, and 
0(*'^(x, n) is the scalar temperature perturbation at position x (n is the direction of light propagation). Using the 
formalism developed in •30] one has 

e(°Hx,n) = J ^E,ef (ryo,k)GO, (24) 
G"e = ^l^Y^oin; ky^^ ^ , (25) 



with respect to a basis where k is fixed. Substituting the above expressions in Eq. ([23]) , and changing basis accordingly 
(c/. HH), one finds 

f d'k (0), , , ,;k.x, 



«^™(x) = J ^er(ryo,k)e'''-r;„(fc; e) . (26) 

Therefore the angular bispectrum is given by (we place the observer in x = 0) 

_ (47^)3(-^)^^+^^+^3 ffkfqfp ^ 

[aiiniiae^ni2ai3,n3) - (2^^ + 1)(2£2 + 1)(24 + I) J (^Tpf ^^imi ejr^^^^ ^9, e^r^g^^l^p, ej {^() 

X (e(°)(ryo,k)e^°>(r;o,q)e(°)(ryo,p)). 
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Remembering Eq. (jlSp and definition (I16p . we see that in order to proceed we need to evaluate the bispectrum of the 
magnetic energy density (pB(k)pB(q)/9B(p)). 



V. THE MAGNETIC ENERGY DENSITY BISPECTRUM 

From the expression of the magnetic energy density given in Eq. pop . we see that its bispectrum is given in terms 
of the six point correlation function of the magnetic field 

1 /■ fj^k d^n d^f) 

W(kW(qW(p)) = (8^y ^^J), ^ (i?.(k)i3.(k- k)i3,(q)i?,(q- q)i?,(p)i3,(p- p)) . (28) 

Since the magnetic field is assumed to be a Gaussian variable, we can use Wick's theorem to decompose the six point 
correlation function into products of the magnetic field power spectrum. To compute the above expression, we then 
use the definition of the magnetic power spectrum Eq. ([T|) and the fact that B*{k) = i?i(— k). Of the total fifteen 
terms obtained using Wick's theorem, seven are irrelevant because they are one-point terms proportional to (5(k), 5(q) 
or (5(p), while eight terms survive. Each of the latter is the product of three power spectra, and consequently contains 
the product of three delta functions (c/. Eq. H])). Two delta functions can be integrated, while the remaining one 
reduces to (5(k + q + p), the homogeneity condition. Starting from Eq. (f28|) . the final result depends on which of the 
variables of the triple integral remains. For example, performing the integration in d^p and d^q and leaving out d^fc, 
one obtains (appropriately renaming the mute indexes) 

(pB(k)pB(q)PB(p)> - Y2^<5(k + p + q) j d3fcP„ (k)F,,(k-k)[i^,i(q + k)+P.i(P + k)], (29) 

where for conciseness we have defined 

Pij{\<L) = Peik) {5,j - kk) (30) 
(note that Pij(k) — _Py (— k)). On the other hand, integrating out d^k and d^p one obtains 

(pB(k)pB(q)ps(p)> = ^^^(k + p + q) j d3qP,^.(q)p^,(q_q)[p^;(k + q) + P,,(p + q)], (31) 

while integrating out d^q and d'^k one obtains 

(pB(k)pB(q)ps(p)) = y^(5(k + p + q) j d3pP^^.(p)p^,(p_p)[p^;(k + p) + P,,(q + p)]. (32) 



This is just a consequence of the fact that the right hand side of Eq. ([28]) is not apparently symmetric under the 
exchange of k, q and p, contrary to the left hand side. Since the final result should be symmetric, we finally set: 

{PB (k)pB (q)pB (P)) = ^^^^^iJi"^^ [j d'~k P^, (k)P,7 (k - k) [Pu (q + k) + Pa (P + k)] 

d^fc (k)P,7 (q - ic) [Pa (k + k) + Pa (p + k)] 



j d3fcp^^.(k)p^.,(p_k)[p^;(q + k)+p^;(k + k)]| . (33) 



Using definitions ([50]) and ^ we have in all generality: 

P,,(k)P,,(q)P,i(p) = A3fc>"g«[(fc • qf + {k ■ pf + [q ■ pf - {k ■ q){k ■ p){q ■ p)] li k < ko , q < kj, , p < k^ . (34) 

and zero else. Due to the complexity of the angular structure and of the integration boundary of the integrals in 
Eq. we cannot derive an exact expression for {pBi}^)pB{^)PB{'P)) which is valid for any configuration of k, q, p. 
We can however give an analytical estimate of the result, which we present in the following. 

We are interested in estimating the behaviour of the integrals in (j33p . From the expression in Eq. [34]it is clear that, 
depending on the value of the spectral index n, the integral could diverge in the infrared limit. On the other hand. 
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Figure 1: The geometrical configuration used to perform tlie integration: k, q and p are free, wliile k is tlie integration 
wave- vector (see p^). 



the angular part always gives a finite contribution. We therefore neglect the angular part for the following estimate, 
and set 



(/3s(k)ps(q)ps(p)) 



(5(k + p + q) 

3847r3 



A' 



d3A;/c"|k-k|' 



|q + kr + |p + k|' 



permutations 



(35) 



To perform the above integration, following [25| . we choose a basis with || k and where the triangle formed by 
k, q, p lies in the plane perpendicular to Cy, in y = 0, see Fig. [TJ We call (p the angle between k and q, cos (j) — k ■ q, 
and a the angle between k and — p, cos(7r — a) = k ■ p. The integration variable k has angles 9 with || k and 4> 
with the plane identified by the triangle formed by k, q, p (c/. Fig. [T|). The angle between k and q is expressed in 
terms of the previously defined ones as 



k ■ q 



sm o cos sm t 



cos cos ( 



(36) 



and the one between k and p is 



k ■ p — — (sin 9 cos 4> sin a + cos 9 cos a) 



(37) 

We remind that the boundaries of the integrals in (j35p are defined by the condition that the momenta coming from 
the power spectrum are bounded by ko- in the first integral for example, the conditions are k < ko, |k — k| < ko, 
|q + k|<fcD. 

Let us first concentrate on the first integral of ([55]) . For negative values of n, it has integrable divergences for k ^ k 
and for k — > — q. We approximate the total result by selecting only these angular configurations, which are the biggest 
contributions to the integral for negative n, and are at least representative of the total result for positive n. By doing 
so, and using the above reference system, the first integral of psp becomes (fi denotes the angular boundary) 



d3fcfc"|k-k|"|q + k|" 



° dk k"+^ I dn 

2n / dkk"+^ 
Jo 



2kk cos 9 



n/2 



q + k- + 2qk{sm 9 cos (j> sin (p + cos 9 cos ( 



|fc - A:|"(g2 + p + 2gA:cos0)"/2 + (/fcZ + p + 2fcfccos(/))"/2|g - k\' 



where in the second equality we have accounted only for the two above mentioned angular configurations: the first 
term of the second equality is the contribution of the angular configuration k ^ k, and therefore 9^0; the second 
one is the contribution of the angular configuration k — > — q, and therefore 9 = it — <j) and (j) — tt. We have inserted 
the factor 27r to simulate the integration in d(j), which should be present at least in the first configuration. We repeat 



n/2 



(38) 
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the same approximation scheme in each term of p5|) . to obtain finally 

( rko p 

\k - fc|"(g2 + p + 2gA:cos0)"/2 + (/fcS + p + 2/cA:cos0)"/2|g _ 




( 



n+2 



n+2 



|fc - A:r(p2 + fc2 - 2pA: cos a)"/^ + (fc^ + P - 2A:A: cos a)"/^ b - fcl" 

|g - ~k\''{p^ + P - 2p~kcos{<j> - + + P - 2gfccos((/) - a))"/2|p _ 



Note that the terms in Eq. ([35]) which share the same wave-vectors collect two by two for the angular configurations 
considered (c./. Eq. (|33p ). It is now possible to evaluate approximatively the above integrals. As already mentioned, 
the apparent divergence for negative n is integrable. Assuming k < q < ko, 'Vfe approximate the first integral in the 
above expression as 



dk k 



n+2 



\k - fc|"(g2 + p + 2gfccos(/))"/2 + (P + P + 2A:A:cos0)"/2|g - k\' 



2 U" /fc" 



dk k 



n+2 



dk k 



2n+2 



dkk 



3n+2 



(40) 



We see that, under this approximation, the angular part plays no longer a role, and the result is the same for the two 
terms of the first line of the above equation. 

Applying the same technique for each integral in Eq. p9p . for the combination /c < g < p < fc^ we find the total 
approximate behaviour: 



(ps(k)ps(q)pB(p)) 



(5(k + p 



(41) 



^ni2n+3 



(n + 3)(2n + 3) 



n j,2n+3 



(n + 3)(2n + 3) 



(n + 3)(2n + 3) 



pnq2n+3 



n 


(3nH 


-3)(2nH 


-3) 




n 




(3n- 


|-3)(27i- 


K3) 




n 




(3n- 


F3)(27i- 


H3) 



3n+3 



-,3ri+3 



3n+3 



7,3n+3 

3n + 3 

i,3n+3 

3n + 3 

j,3n+3 

3n + 3 



for k < q < p < ki 



while if g < A; < p we have to exchange k and q in the above expression, and so on with all the ordered permutations 
of the wave-numbers. 

Observing Eq. (j4ip . we can confirm what pointed out in [25], i.e. that there are two distinctive spectral regimes for 
the bispectrum. For fiat and blue magnetic field spectra, with n > — 1, the infrared limit fc of the bispectrum is 
white noise; Eq. (|4T|) is in fact dominated by the constant terms k^~^^ / {3n + 3) . On the other hand, for red magnetic 
field spectra n < —1, the bispectrum is divergent in the infrared limit. As we will see in the next sections, the 
divergence can go as fc^"^'^ or as k^'^'^^, depending on the wave- vector configuration. The same behaviour holds for 
the magnetic energy density power spectrum, but in this case the discriminating value is tz = —3/2, and the infrared 
divergence for n < —3/2 goes as fc^""*"^ [13, [2§| (in [23|, it has been found that for the limiting value n = —3/2 the 
white noise spectrum acquires a logarithmic dependence on k: this is the case also here for the corresponding limiting 
value n = —1). The above approximated result is valid only for k, q and p smaller than the magnetic upper cutoff 
ko, while in general they do not need to satisfy this bound. As already mentioned in section [III in [131 it has been 
found that the magnetic energy spectrum goes to zero at fc = 2kD, due to the convolution boundaries (c./ Eq. (jl2p ). 
As we will see in the next sections, the same behaviour holds also for the bispectrum (this is verified exactly in the 
collinear configuration) . 

The above equation (j4ip is a general approximation to the magnetic field energy density bispectrum in the infrared 
limit. We now compare it with the result coming from a specific configuration of the wave-vectors, the collinear 
configuration, for which we have an exact result. We find that the above expression can be considered quite a good 
approximation to the true magnetic field bispectrum in the infrared limit. We also give explicit formulas for the 
squeezed and equilateral configurations, for which, however, we do not calculate the exact result. We find that 
the three configurations give a comparable white noise contribution for n > — 1, while if n < —1 the collinear and 
equilateral configurations diverge in the infrared limit as k^"'^^, while the squeezed one diverges as /c^""*"^. 
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A. Collinear configuration 



The collinear (or flattened) configuration is given by two equal wave-vectors, while the third one points in the 
opposite direction: for example, p = q and k = — 2q. In this case, it is possible to calculate the bispectrum ([55)1 
exactly. The three permutations of k, q, p of this configuration should be present in the symmetric expression (j33l) : 
this gives in the end 



(ps(k)pB(q)ps(p)) 



collinear 



Puik + k) + Pa ( --k 



+ P,,(k-k)P,, --k 



(42) 



Therefore, for the collinear case we find the following expression, using Eq. ([M)) : 
(ps(k)pi3(q)PB(p))| 

collinear 



(5(k + p + qy 3 







n 




[ (fk k" 






k + k 


Ivi 









, A" X 
5767r3 

(fc-k + 2fc)2 (A:.k + fc)2 



(fc2+3k-k+2fc2)(fc2 _ (fc.k)2) 



k + k 



k + k 







n 












IV2 









(fc-k + 2fc)2 (fc.k-2fc)2 (fc2-4fc2)(A:2_(/t.k)2) 



l+k 



is - k 

2 ^ 



16 



l+k 



is - k 

2 ^ 



+ k^p + k^q}, 
where Vi denotes the volume given by the three conditions 



(43) 



and V2 is given by the conditions 





k 


< 


ko 


|k/2 


-l-kl 


< 


ku 


|k 




< 


ko 




k 


< 


ko 


|k/2 


+ k| 


< 


ko 


|k/2 


-k| 


< 


ku 



(44) 



(45) 

The last term of Eq. 221 PijO^)PjiXk — k)Pii — k^ , becomes equal to the first one by changing k to — k and the 
integration volume accordingly. 

It is possible to calculate Eq. ([43]) exactly for the selected values of the spectral index n = 2 and n = —2. This 
is due to the fact that, in this configuration, the integration over the angle (j) becomes trivial (c./. Fig. [T]): since 
p = q and k = — 2q, the integrands in Eq. P5)) depend only on cos^ and the boundaries given by Vi and V2 can be 
made explicit with little difficulty. The details of the calculation are given in Appendix [Bl while the result is shown 
in Fig. [21 In the case n = —2 the calculation is quite involved, therefore we have evaluated only the infrared part, up 
to fc < fc£)/2. On the other hand, the case n = 2 is simpler, and in this case we found a general, exact expression. 
This expression confirms that the cutoff of the bispectrum is at k = 2kD, as we would expect from the analysis of 
the spectrum (see Fig. [3]), and as can be viewed easily from the last inequality of the boundary conditions of Vi, 
|k + k| < kr), which shows that the maximal allowed value for k is 2k d. 

Knowing the exact result, we can test the goodness of the approximation given in the last section at least in this 
configuration. Reducing the general result of Eq. ([1T|) in the collinear configuration, we find: 



{pB(k)pB{q)pB{p)) 



collinear 



(5(k + p + q) 
1447r2 

+k^p+k 



A' 



23"+3(2n + 3) 



2"+i + 1 2 



3n+4 



1 



3n + 3 



9} 



(46) 



Given that k < kjy, if n < —1, this expression is divergent for fc ^ as fc3"+3^ while for n > —1, is it white noise. 
Consequently, the case n — —2 exhibits a divergent behaviour as k~^ , while the case n = 2 is regular, as can be seen 
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n=-2 

n=2 




Figure 2: The magnetic field bispectrum in the coUinear configuration p = q = — k/2, normahsed by the quantity 
A^k^^"^ /{576tv^), as a function of k/ko, for n — 2 (left plot) and n — —2 (right plot). We only show the infrared region 
k < kol'i.. The blue, solid line is the exact result, while the red, dashed line the approximation given in Eq. (|46p . 



n=2 n=-2 




Figure 3: The magnetic field bispectrum in the coUinear configuration p = q = —k/2 (blue, solid) and the magnetic field 
spectrum to the 3/2 (red, dashed), both multiplied by the phase space density k"^, as a function of k/kn for n = 2 and n = —2. 
Note that in the n = —2 case, we only calculated the bispectrum up to fc = kD/2, while the spectrum is known up to = ko- 

in Fig. [21 where the true and approximated result are compared. In the regular n = 2 case, the bispectrum is not 
pure white noise but shows a mild dependence on k: our approximation does not capture this dependence, but only 
the infrared white noise behaviour. In both cases, our approximation underestimates the true result by a factor of 
two. 

In Fig. [21 we compare the exact result of the bispectrum in the collinear configuration with the magnetic spectrum 
to the power 3/2, for n — 2 and n = —2, both multiplied by the phase space density {k/koY ■ For n = 2, they are of 
the same order of magnitude, as one would expect. For n — —2, the bispectrum goes as fc'^""'"'^, while the spectrum 
as fc^"+'^. The spectrum approaches the bispectrum amplitude as k grows, however, the exact bispectrum has been 
calculated only for k < A;d/2, and this is the region shown in the plot. 

B. Squeezed configuration 

In the squeezed configuration one wave-vector goes to zero while the other two are equal but opposite in direction. 
Expliciting the case in which q ~ 0, k = — p, from Eq. ([55]) we have 

+ Pji (q - k) [Pa (k - k) + Pu (k + k)] 
+ Pji (k + ic) [Pu (k + k) + Pu (q + k)] 

+ (q^p~0,k^q) + (q^k~0,k^p)} . (47) 
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Using Eq. I34[ and grouping the terms which are mutually equal, we find: 
{pB{k)pB{q)pB{p))l 

QUCCZ( 



,5(k + p + q)^3^ 



3847r3 



(48) 



2 

3< 







n 






k-k 




q + k 


'Vi 









(A;-k-fc)2 {k-ti + k)^ (k-q-k-k + q-k-P)[k-q- (fc-k)(fc-q)] 





2 ' 






k k 






hk 



k k 



q + k 



d^~k ~k'' 



V2 



d-^k A:" 



k k 



q-k 



2n 



{k-k-ky (fc2_2k-fc + fc2)[fc2-(/j.k) 



k k 



k-k 



k-k 



(fc-k-fc)2 {k-q-kf (k-q-k-k-q-k + fc2)[k.q- (fc.k)(fc-q)] 



'^3 k k 

+ (q^p~0,k^q) + (q^k~0,k^p)} 
where Vi is given by the conditions 



q k 



+ 



k-k 



q-k 



V2 by the conditions 



and V3 by the conditions 





k 


< 


ko 


k- 


k| 


< 


ko 


q + 


k| 


< 


ko 




k 


< 


ko 


k- 


k| 


< 


ko 




k 


< 


ko 


q- 


k| 


< 


ko 


k- 


k| 


< 


ko 



(49) 



(50) 



(51) 

We do not have an exact calculation of the bispectrum in the squeezed configuration. This is due to the fact that, 
contrary to the coUinear case, the integration over the angle (j) is not trivial (c./. Fig. [T]). For example, for the 
case q ~ 0, k = — p, the angle (f> — > tt/2: therefore, taking for example the first integral in Eq. (|48|) . we see that 
it contains the term {q^ + fc^ + 2kq sin 9 cos <]))"' ^"^ , and the integration boundary over 9 and (j) given by Vi becomes 
very complicated. Having already an estimation of the goodness of our approximation in the coUinear case, we do 
not dwell on the calculation for the squeezed configuration, and use only the approximated formula. Observing 
the boundary conditions given by the integration volumes Vi, V2 and V3 we can however confirm that, also in the 
squeezed configuration, the bispectrum goes to zero a,t k = q ^ 2ku, since these are the maximally allowed values for 
the wave-numbers. 

We reduce the general expression given in Eq. ([1T|) in the squeezed configuration, and find: 



(PB(k)pi3(q)ps(p))ls 



5(k + p + q)^3 



1447r2 



2n 



(n + 3)(2n + 3) 



6n{n + 2) 



-k 



3n+3 



71 + 1 



+ (g ^ p ~ , fc ^ g) + (g 



(3n + 3)(2n + 3)(n + 3) 
k~0,k^p)\ . (52) 



For n > — 1, the resulting white noise plateau has the same amplitude as in the coUinear case. However, for n < — 1 
the divergence for g ^ is g-^n+s. therefore, it is weaker than in the coUinear case, and reaches the coUinear behaviour 



.,3ri+3 



only in the limit k ^ q ^ 0. 
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C. Equilateral configuration 



In the equilateral configuration the wave-vectors form an equilateral triangle. With q — kq, and p = kp, using 
Eq. (j33p and regrouping the equal terms one gets 



{pB{k)pB{q)PB{p))\ 



equilateral = ^^^^^^^F^ ^ / d'kP,,{k){p,,{k-k)[PAkq + k)+PAkp + k)] (53) 

+ P-ji {kq - k)Pii {kp + k) + (k ^ q,kq ^ qk ,kp ^ qpj + (k ^ p , kq ^ pq , kp ^ pkj | . 



Using Eq. (|34p . we can rewrite the above expression explicitly as 



{pB{k)pB{q)pB{p))\ 



^(k + p + q) ^3 



equilateral 



3847r3 



X 



(54) 



d^k /c" 



Vi 



k-k 



d^k A:" 



Vi 



d-^k ¥• 



V3 



k-k 



kq ~k 



kq + k 



kp + k 



kp + k 



(fc-k-fc)2 {k-kq + kf {k-kq-k-k + kq-k^P)[k-kq-{k-k)(k-kq)] 



k k 



kq + k 



+ 



k-k 



kq + k 



(fc-k-fc)2 {k-kp + kf {k-kp-k-k + kp-k-k'^)[k-kp-{k-k){k-kp)] 



k k 



kp + k 



2 _ 2 

k — k kp + k 



{k ■ kq — kY (k ■ kp + fc)^ {kq ■ kp — kq ■ k + kp ■ k — k^)[kq ■ kp — {k ■ kq){k ■ kp)] 



kq — k 

[k ^ q,kq ^ qk ,kp ^ qpj + [k^ p,kq ^ pq 
where again Vi is given by the conditions 



kp + k 
, kp -> pkj I 



kq — k 



kp + k 



(55) 



k < ki) 
\k-k\ < ko 
\kq + k\ < ko, 



(56) 



and similarly for V2 and V3. In this case as well, we cannot solve the above integrals exactly. Like in the squeezed 
configuration, the integration in d(j) is non-trivial, since (j) — 27r/3 and, for example, the first integral of Eq. 
contains terms like [k"^ + fc^ 
Eq. (|4ip . which gives simply: 



2fcfc(^ sin 9 cos (f) — \ cos 

6'))"/2. We therefore use the approximated expression in 



(ps(k)/3i3(q)pi3(p)) 



<5(k- 



q) 



equilateral 



1447r^ 



A' 



6n 



u3n+3 



(n + 3)(3n + 3) 



i,3n+3 

n + 1 



+ ik^q) + ik^p)\ . (57) 



For n > —1, we find again a white noise plateau of the same amplitude as in the other configurations; for n < —1 
the divergence for fc is the same as in the collinear case. At first sight, this result might not seem correct: in 
the collinear case, in fact, by definition the wave- vectors are collinear and therefore the limits k ^ k and k — > — q 
collapse into a single wave-vector configuration. In the equilateral case, on the other hand, they do not: we would 
therefore naively expect the same infrared behaviour of the squeezed configuration. However, the infrared divergence 
occurs for k — q ~ p 0, and in this limit k — > k and k — > — q are no longer distinct. Therefore, we do expect a 
^3n-i-3 behaviour also in the equilateral case, equivalent to what we find in the collinear case and also in the squeezed 
one when we let not only g, but also k ^ (c ./ Eq. 

We can conclude that, although it neglects the angles, the approximation in Eq. (|4T|) does recover the correct 
behaviour of the bispectrum in the analysed configurations. However, neglecting the angles certainly introduces an 
inaccuracy, because one does not account precisely for the weight with which the different configurations contribute 
to the total result. We were able to compare the approximated result with the exact one only in the collinear 
configuration, and we found an underestimation of a factor of two both for negative and positive spectral indexes. 
However, this does not ensure that the total, exact bispectrum is altogether only a factor of two higher than what 
given in Eq. (j4ip . neither that it has exactly the same dependence on wave-numbers when we significantly deviate 
from the infrared limit. 
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VI. THE CMB BISPECTRUM 



Given the magnetic energy density bispectrum ((0B(k)/9B(q)/3B(p)), we can now evaluate the CMB bispectrum 
Eq. (I??)) . We use the approximated magnetic energy bispectrum Eq. (|¥T|l . 

^3t,3«+3 

(pB(k)ps(qW(p)> ^ 6ik + p + q)—J^J{K,Q,P) (58) 
^ ' (n + 3)(2n + 3)^ (3n + 3)(2n + 3) ^ (n + 3)(2n + 3) 

J p3n+3 J pnQ2n+3 ^ p3ri+3 



(3?i + 3)(2n + 3) (n + 3)(2n + 3) ^ (3?i + 3)(2n + 3) 

+ for K <Q <P <1, 

n + 1 

where K = k/ku and so on denote normahsed wave-numbers. We want to estimate the reduced bispectrum bi-^i^i^ 
introduced in [421 



)=C™&W3, (59) 

where Gg^l^^"^^ is the Gaimt integral. We use the procedure described in [i^: starting from Eq. ((27)l . substituting 
in it Eq. ITSl) and Eq. (fT6|) . and using expression ((58| for the source, we find: 

Tra^yl^A;^"^^ f°° 

bi,i,i, = . / / d/Cif^ / g2 / dP P^je,{Ky)je,{Kx)n,{Qy) 

O/'rol "'0 "'0 Jo Jo 

J,, (0x)j,3 (Py) (Px) X(A', g, P) , (60) 

where y = fcuryo and a; = /cd'', and r comes from the decomposition of the delta function in ([55]) (see [l^). Using the 
definition of the bispectrum T{K, Q, P) given in (|58)) . the above equation becomes 



/ dxx^ dKK'u,iKy\u,iKx) / dgQ^^,, (Q2;)j,, (Qx) / dPP^ j,,{Py)j,,{Px) 
Jo Jq Jo Jo 

X {a{n) [i^"Q2n+3 _^ ^np2n+3 _^ Q«p2n+3j _^ ^(^^^ [2ii:3"+3 + g3«+3] + c(7l)} 

+ permutations , (61) 

where a{n) — n/{n + 3)/(2n + 3), b{n) — n/{3n + 3)/(2n + 3), c{n) = l/{n + 1), and one adds the six ordered 
permutations of Q and P which entail permutations of £i, £2, is- 

In order to estimate the bispectrum, we substitute the upper boundaries in Eq. (|6ip with the interval [0, 1] in all 
the integrals over the momenta, since the Bessel functions peak at very low momentum: je^{Py) peaks at P ~ ^3/2/, 
and y 3> 1. Because of the form of the source T{K,Q, P), in ([6T|l at least one integral over the momentum is not 
influenced by the source. Following [44], for each of these integral we use the approximation (c/. Eq. 6.512 of [ist ) 



dP P^u, {Py)3,, [Px) ^ \ ^^^^ , (62) 



we then solve the integral in dx using the delta function and obtain for the first term for example, 

a{n) 



f dK K^+^jI (Ky) f dQ Q'^+'jI (Qy) ., 
Jo Jo 



(63) 



and so on. Approximate expressions for this kind of integrals are discussed in appendix 1X1 

If n > — 1, in Eq. (|6ip we retain only the white noise term c{n). All permutations give the same result in this case, 
and we find finally 



7r7a3(n + 3)3 (P2)3 1 

Ki,t, - Qfi ' 1 3 71 ^' forn>-l. 64 

96 n + l p;*^! (korio) 

Values of the spectral index n < — 1, for which the source is not pure white noise, are a bit more involved. As in 
the spectrum case (c/. section HIH) . we cannot give a general expression valid for every n < — 1, since the way to 
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approximate integrals like those in ((63)) depends on the actual value of the power law exponent. Therefore, we give 
explicit expressions only for two values of the spectral index: n = —2, and n — > —3. Fixing the spectral index to 
n = —2, one finds 



288 {komfXh 
permutations , 



log 



for n 



2fcz???o 1 
3tt £i 



1 

T2 



1 , / kom 



kDVo 1 



37r 



(65) 



It is important to remark that the squeezed limit of the above expression must be taken with £3 <^ £2 — ii, since this 
expression has been derived from the wave-number configuration P < Q < K . We see that in this case, the dominant 
term in the bispectrum (of the order log (fc^ 770/^^3)) correctly corresponds to the one coming from the dominant term 
in wave-number space, p2n+3 ^pj^^, permutations must be treated accordingly: for example, for Q < P < K one has 
£2 « £1 ^ £3. 

For n —3, we solve the integrals setting n = —3, therefore using approximation l|A3|) with m = —1, and we find 
then 



TT^ n(n + 3)2 (^2)3 



288 2n + 3 
permutations , 



1 

W2 



1 

for n ; 



1 

-2*3 
s -3. 



TT n + 3 
16n + 1 



koVo 



1 



1 



(66) 



where the same considerations as above apply for the squeezed limit. The second term in the above expression, 
coming from the term proportional to h{n) in Eq. ([61]) is sub-leading, since it contains a factor n -|- 3. Note that since 
{B'^) oc (n + 3)^^, the leading term of the bispectrum diverges in n ^ —3 as [n + i]-^, like the spectrum (c/. Eq. ((22l) l: 
this divergence is connected to the infrared divergence of the magnetic energy ^ . 

The leading term of the above result Eq. ((66|) . reduced to the squeezed and the equilateral configurations, gives the 
same result as found in [28] (c/. eqs. (17) and (18) and discussion thereafter, we remind that we use a — 0.1). 



VII. ESTIMATION OF THE SIGNAL 

Since the signal-to-noise ratios {S/N) we will be interested in is some function of the maximum multipole a given 
experiment can reach, ^max 3> 1, we can use the flat-sky approximation (46l. |47| and write for the bispectrum 



(a(fi)a(f2)a(4)) = [2nf5'^''\£i23)B{£i,£2,£3), 



(67) 



where ^123 = ^1+^2+^3- With this notation, the reduced bispectrum hg^g^g.^ coincides with the bispectrum B{£i,£2, £3). 

Our goal now is to quantify the level of NG coming from the stochastic magnetic field and eventually to give a 
bound on the amplitude of the magnetic field. One way to do it is to exploit the present bound on the primordial 
local non-Gaussianity parametrized by the quantity /^l- we mentioned in the introduction, the search for a 
non-vanishing bispectrum of a local type has given so far a null result and currently is bounded in the range 
—9 < < 111. As the shape of the non-Gaussian signature from the stochastic magnetic field may be different 
from the one of the local type, one may not directly apply the bounds coming from WMAP5 whose search for non- 
Gaussianity is optimised to search for local primordial contribution. Instead, we proceed in the following way. First, 
we define the Fisher matrix (see, for example, (i^ l 



/sky 

(27r)27r 



(f£i(f£2(f£3S^^Hh 



23) 



B^{£i,£2,£3)Bj{£,,£2j3) 
6C{£i)C{£2)C{£3) 



(68) 



where /sky is the portion of the observed-sky in a given experiment and i (or j)— (mag, loc). The first entry i^mag,mag of 
the Fisher matrix corresponds to the signal-to-noise ratio {S/N)'^ provided by the stochastic magnetic field to the non- 
Gaussianity. We have defined the power spectrum in the flat-sky approximation by (a{li)a{l2)) = (27r)^(5(^) (/i2)C(£i) 



^ The above expression is valid only for n —3 so the denominator is always finite. Note however that the apparent divergence for 
n = —3/2 is just an artefact due to our approximation (c/, Eq. ((TSJ): n = —3/2 would correspond to a threshold value for which 



|P-B(fc)P diverges logarithmically for k ^ and is not simply white noise. 
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with £^ C{t) = A/n and A ~ 17.46 x 10~^ is the ampUtude of the primordial gravitational potential power spectrum 
computed at first-order. In other words, we assume that the two-point correlation function is dominated by the usual 
adiabatic contribution from inflation. Finally, the local bispectrum is given by [iHl 

o ^loc /la / 1 \ 

i3ioc(^i,^2,4) = -^^^ [jr(2 + cyci.j . (69) 

Notice that all these expressions are obtained in the Sachs- Wolfe approximation. We will return back to this point 
shortly. 

Next, we define an effective /^^ which minimises the defined as 

r ( /§i-Bloc(^l,^2,^3)Lloc^l ~ ^mag(^l,4,4) 

X jat,ai2ai30 (I123) 6C(£i) ^(^2) 



One finds 



reS _ ^mag.loc 
ioc.ioc 



(70) 



The signal-to-noise ratio for the primordial local case has already been computed in the flat-sky approximation in 
Ref. [46]. The result is that -Flocjoc - {4=/ T^^)f sky A{f ^l)'^ £l,^^ log(Cax/Cin)- The logarithm is typical of scale 
invariant power spectra and ^^in is the minimum multipole compatible with the flat-sky approximation. The physical 
meaning of /§l is the following: it is the best value of the local f^f^ which best mimics the bispectrum from a 
stochastic magnetic field background. As such, we can apply to this value the current observational limits. 

We start with the simplest case n « —3. Indeed, for n close to —3, the leading term of the bispectrum is of the 
same form of the local primordial bispectrum ([M)) in the squeezed limit £3 <C -^i — £2 and we immediately find 

^^^" 288^2 2n + 3 p^, " ^" + (, (lO-^Gauss)^ j ' forn « -3 . (71) 

In the case in which the bispectrum is independent from the multipolcs, that is for ?i > — 1, we find 

.^a^ (n + 3)3(i?^)V wy 1 -^6^,10-7(^ + 3)^ ( (^^) V for , > _i 

•^NL- 2304^2 „+i ^3^^ ^g^J log(Cax/Cin) " n + 1 ^ (10-9Gauss)2 ^ ' ^ ■ 

(72) 

Finally, for the case n = — 2, we find 

^^1^- 2304^2^?,, [£d) log( W^min) - ^ " (lO-^Gauss)^ j ' ^""^ " " (^"^^ 

In all numerical estimates we have taken £0 — kor]o — 3000, ^max ~ 750, ^min ~ 10, a ~ 0.1, and Eq. (El). We see that 
the effective value of non-Gaussianity is smaller than the present upper bound of 0(10^) on f^l^ \2\ for magnetic 
fields 0(10) • 10"^ Gauss for n « -3 and 0(20) • 10"^ Gauss for the other cases^ 

Accounting more precisely for the value of the damping scale fc^) as a function of the spectral index and of the 
magnetic field amplitude using Eq. ([7]), we obtain 



B^) < 9 nGauss for n = -2.9 



'B^) < 25 nGauss for n = -2 



< 20 nGauss for n = 2 (74) 



We have obtained similar estimates repeating the same procedure to define an efli'ective non-Gaussianity parameter starting from a 
primordial equilateral configuration for which WMAP5 limits exist. In such a case the primordial equilateral configuration is peaked 
for £1 £2 ^ £3 and the effective non-Gaussianity parameter scales with £max with one power less than the corresponding one obtained 
from a local primordial bispectrum. 
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The corresponding bound on the magnetic field amplitude B\ {cf. Eq. ([5])) on the scale A = 0.1 Mpc is unchanged 
for n —* —3, it becomes Bx < 26 nGauss for n = —2, and is less stringent as n grows, becoming irrelevant for n = 2: 
B\ < 2 /iGauss. This is a consequence of the fact that the procedure of using an effective /nl returns a bound on the 
integrated magnetic field spectrum, and therefore for very blue spectra the constraint on large scales is irrelevant. 

A word of caution is in order here though. In all our estimates, we have used the Sachs- Wolfe approximation 
for all bispectra. This is certainly a sufficiently good approximation for an experiment like WMAP whose maximum 
multipole is £niax ^ 750. This is because the transfer functions for both the scalar contribution to the CMB anisotropics 
from the stochastic magnetic field and the one from the inflationary adiabatic modes may be taken roughly equal to 
unity up to £ ~ 750 and they do not affect the computation of the Fisher matrix elements, see [4^ [4^. However, for 
higher multipoles, say i ^ 2000, typical of an experiment like Planck, the inclusion of the transfer functions will be 
crucial because the anisotropics from the adiabatic inflationary modes get an exponential suppression due to the Silk 
dairming, while the ones from the scalar modes from the stochastic magnetic field show a much milder suppression 
[27II4O.]. This will increase the value of /^l. Needless to say, the inclusion of the vector and tensor contributions from 
the magnetic field will help to increase the non-Gaussian signal too. 

While writing this paper, the preprint [2^ appeared where the computation of the bispectrum from a stochastic 
magnetic field background was presented for the case n « —3. Our findings agree with those in Ref. [28j and extend 
them to other values of the spectral index and by the estimation of the signal-to-noise ratio and of the effective 
non-Gaussianity parameter. 
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Appendix A: INTEGRALS OF BESSEL FUNCTIONS 



In order to evaluate both the magnetic field spectrum and bispectrum at large angular scales, we need to evaluate 
integrals of the type 



(Al) 



with y ^ 1. This integral can be expressed generically in terms of hypergeometric functions; however, good approxi- 
mations can be found, which are much simpler. 

For TO = 2, the integral can be performed exactly: one has 



(A2) 



where since y ^ £ we used the expansion of the Bessel functions for large arguments. 

For TO < 1, the integral reaches a constant value for y ^ £, and can therefore be evaluated in the limit y ^ 00. We 
find 



dxx^^jfix) 



1 r v^r(i^)r(£ + ^) 
4 [ r(i-f)r(£ + 2^) 

for TO < 1 , y ^ i 



Jy_ 

TO — 1 



-I- sin(7r^ - 2y) 



4r(i - f ) 



(A3) 



The case m = 1 is a bit more involved: the integral (|A1[1 grows logarithmically with y and cannot be evaluated 
with the same approximation as before. In this case we set 



dx / 7r 7r \ 1 
dxxjjix) - — cos^ i^x~-£--j ~- [log(y) - log(£)] for y > . 



(A4) 



We are neglecting the subdominant contribution to the integral of the interval [0,£], therefore this approximation is 
slightly underestimating the true result. However, it captures the correct behaviour in £ and y. These approximations 
are shown in Fig. [4l 
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m=2 m < 1 




Figure 4: The approximations for the integral in Eq. (IA1[I . Upper left plot, for m — 2: the integral (solid) and the approximation 
y/2 (dashed) are shown for £ = 20, £ — 100, £ — 500 as a function of y. Upper right plot, for m < 1: the approximations for 
y ^ £ (solid) and for £ 1 (dashed) given in Eq. (IA3P are shown as a function of £ for m = 0, m = —1 and m = —2. Lower 
plots, for m = 1: the integral (solid) and the approximation in Eq. (|A4|) (dashed) are shown as a function of j/ for £ — 20 (left 
plot) and as a function of £ for y = 100 (right plot). 



Appendix B: BISPECTRUM IN COLLINEAR CONFIGURATION 

In the following appendix we describe the technique used to calculate the magnetic energy density bispectrum in 
the collinear configuration Eq. (|42p . 

Due to the complexity of the calculations we restrict to analytical solutions of the bispectrum integral for two 
representative spectral indexes: the case n — 2 (the typical spectrum of a magnetic field generated by a causal 
mechanism), and the case n — —2 (in order to investigate the behaviour of the spectrum also for negative spectral 
indexes). 

The magnetic energy density bispectrum in the collinear configuration is given in Eq. (j42p . From this we extract the 
integral in the momenta which is given by the permutation over the three momenta K, P, Q of three basic integrals^: 

I{K)= f dK f dx{Ia{K,K)+h{K,K) + UK,K)), (Bl) 



For simplicity of notation in this appendix we use re-scaled variables: K = k/k]j, Q = q/k]j, P = p/kjj and K = k/kjj 
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where x = k-k. The functions Ia{K, K), h{K, K), Ic{K, K) are: 

K'^ - ~ -n. 

Ia{K,K) = K^+'^i— + + KKx)^{K^ + + 2KKx)-^+^ 



/Sk'^ + 24:Kk^x + K^{1 + x^) + 3K^kx{3 + x^) + K'^k'^{7 + I9x^)\ 
V K"^ + 4^2 + AKkx / 



ib{K,k) = k'^+'^i^ + k'^ + Kkx)''/'^{k^ + ^K{K-Akx)Y''^ 

/(32if4 + AK^k'^il - 5x2) + ^4(1 + a;2))\ 



((^2 +4i^2)2 _ 16i^2X22.2) 

^2 + i^2 _ 2Kkx)-^+'i {k^ H 
(8^:^^ - 24Kk^x + K'^il + x^) - 3K^kx{3 + x^) + K'^k'^{l + 19a;2))^ 



Ic{K,k) = k^+'^iK^ +k^ -2Kkx)-^+^{k^ + ^K{K -Akx))""^^ 



(^2 + _ iKKx) 



(B3) 



(B4) 



Wc note that due to the symmetry k — !■ —K we have that the first and the third integrals are indeed the same 
Ia{K) — Ic(K), therefore to obtain the energy density bispectrum in the colhnear configuration we need to solve only 
the two integrals of Ia{K, k) and Ib{K, k). 



1. Integration Domains 

The sharp cut-off of the PMF spectrum at the damping scale kn, imposed to account for the magnetic fields 

suppression on small scales, leads to many conditions on the angle K ■ K. This causes the integration domain to be 
split into various sub-domains. The conditions are different for /„ and /b, therefore for simplicity in the following we 
consider the two integrations separately. 

2. Domains of la 

The sharp cut off imposes: 

k < 1 

i— + + KKx) < 1 
{K^ +k^ + 2Kkx) < 1 
This leads to the following integration scheme: 

1) Q<K <l 

n\ — K p\ p\ p— — — — ^^-^ — 

/ dk dxia{k,K)+ dk dxIa{k,K) 

Jo J-1 Jl-K J-1 

2) 1<K <2 

r - r 

I dK dxIa{K,K) (B5) 

Jk-1 J-1 
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3. Domains of 

The sharp cut off hiiposes: 

k < 1 

(— + KKx) < 1 

K^ ~ 

(—+K^ - KKx) < 1 
This leads to the foUowing integration scheme for < if < 2: 



f ' dk f dxh{k,K)+ I ' dk I _dxh{k,K) 

Jo J-1 J-1 + k2^4 + K± 

(B6) 



in the interval ^ < K < 1 the integral collapses to zero. 

4. n=2 

First we consider the case n — 2 which is the easiest from the point of view of the calculations. In fact the angular 
integrand functions for this spectral index simply reduce to: 

Ia{K,k,x) = ^k^{8k^ + 24Kk^x + K^{l + x^) + 3K^kx{3 + x^)+K^k^{7+19x^)) 

Ib(K,k,x) = —k'^{32k^ + 4:K^k^(l-5x^)+K'^(l + x^)) (B7) 
16 

Once performed the angular integrations, following the integration scheme reported in the previous paragraph, the 
radial integrations become trivial and the result is: 

^(^)l'-^^(3-3^+-r-^ + — -l6 + 256 - 5376o) ^^^^ 

In Fig. [2] we have shown the result for n = 2. We note that, as it happens for the energy density spectrum, also the 
PMF energy density bispectrum goes to zero for if = 2 as expected. 



Here we consider the case n = — 2. The functions la and for this spectral index reduce to: 

4(8A"* + 24:Kk^x + K^{1 + x^) + 3K^kx{3 + x^) + K^k^{7 + 19x^)) 



Ia{K,K,x) 



{K'^ + K^ + 2KKx)^{K'^ + AK"^ + AKKx) 



2 



^ m{?,2K^ + AK-K\l-^x^)+K\l + x^)) 

We note how these functions are far more complicated than the ones for the n — 2 case. Once performed the angular 
integrations in both the integrals we have the appearance of absolute values like \K — 2K\ and |if — -ftTj, their presence 
influences the integration domains creating further splitting into several sub-domains. Since we are interested in the 
effect on CMB where only the low K part of the spectrum has a role we restrict ourselves to the K < 1/2 region of 
the spectrum. The analytical result for n — —2 unfortunately has a very long and complicated form, therefore, for 
the sake of simplicity, we show only the infrared limit: 

24.674 



UK) 
h{K) 

l{K)r. 



X3 

24.674 

73.8367 
K^ 
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Fig. [2] shows the exact result. 
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